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Abstract 

Using an approach based, amongst other things, on Proposition 1 of Kaluza (1928), Goldie 
(1967) and, using a different approach based especially on zeros of polynomials, Steutel (1967) 
have proved that each nondegenerate distribution function (d.f.) F (on R, the real line), 
satisfying F(0— ) = and F(x) = F(0) + (l-F(0))G(x), x > 0, where G is the d.f. corresponding 
to a mixture of exponential distributions, is infinitely divisible. Indeed, Proposition 1 of Kaluza 
(1928) implies that any nondegenerate discrete probability distribution {p x : x = 0, 1, . . .} that is 
log-convex or, in particular, completely monotone, is compound geometric, and, hence, infinitely 
divisible. Steutel (1970), Shanbhag & Sreehari (1977) and Steutel & van Harn (2004, Chapter 
VI) have given certain extensions or variations of one or more of these results. Following a 
modified version of the C.R. Rao et al. (2009, Section 4) approach based on the Wiener- 
Hopf factorization, we establish some further results of significance to the literature on infinite 
divisibility. 
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1 Introduction 



Several papers, including those of Goldie (1967) and Steutel (1967), address the preservation of 
infinite divisibility property under mixing. Steutel (1970) and Steutel & van Harn (2004) have 
reviewed and unified most of the literature in this respect. This consists of, amongst other things, 
the result, due to both Goldie (1967) and Steutel (1967), that each nondegenerate d.f. F (on R, 
the real line), satisfying F(O-) = and F(x) = F(Q) + (1 - F(Q))G(x), x > 0, where G is the 
d.f. corresponding to a mixture of exponential distributions (or, in other words, corresponding 
to a density with its restriction to (0,oo) completely monotone), is infinitely divisible (i.d.), and 
some of its extensions and variations. It may be worth pointing out here that the result relative to 
exponential distributions mentioned above is usually referred to as the Goldie-Steutel theorem or 
Goldie-Steutel result; to be short, we refer henceforth to this simply as GSR. 

Proposition 1 of Kaluza (1928), a renewal theoretic result, implies clearly that any log- 
convex or, in particular, completely monotone, nondegenerate discrete probability distribution 
{p x : x = 0, 1, . . .} is compound geometric, and, hence, i.d. In view of the information in Steutel 
(1970, p. 89) (and the closure property of the class of i.d. distributions), one can then see that, 
as a corollary to Kaluza's result, any nondegenerate d.f., F (concentrated) on [0, oo), that is 
differ entiable, with log-convex derivative, on (0, oo), is i.d. (In Section 2 to follow, we are to 
discuss some basic properties of log-convex sequences and log-convex functions, of relevance to this 
study.) The two results mentioned here obviously extend GSR and its discrete version (respectively, 
in the reverse order). In his argument to obtain GSR, Goldie (1967) also uses the proposition referred 
to, but, he does so in conjunction with a certain characteristic property of the i.d. distributions 
on [0, oo), based, essentially, on mixtures of Poisson distributions. (We view here the degenerate 
distribution at the origin as Poisson with mean zero.) Incidentally, the characterization, just 
referred to, shows especially, that GSR implies its discrete version and vice versa (and that this is 
also so for Theorems VI. 3. 10 and VI. 7.9 of Steutel & van Harn (2004), respectively). 

Shanbhag & Sreehari (1977, Theorem 5) and Steutel & van Harn (2004, Theorems VI.3.13 and 
VI. 7.10) provide us with extended versions and variations of GSR; in these contributions, Theorem 
2.3.1 of Steutel (1970) has implicitly played an important role or provided a motivation. Also, it 
may be noted here that Theorem VI.3.13 of Steutel & van Harn (2004) is a corollary to Theorem 5 of 
Shanbhag & Sreehari (1977), and that some applications of GSR and its extensions have appeared 
especially in Shanbhag & Sreehari (1977) and Shanbhag et al. (1977). The main theme of the 
arguments used by Steutel (1967, 1970) and Steutel & van Harn (2004) to obtain the results implied 
is obviously based on zeros of polynomials, but, more recently, C.R. Rao et al. (2009, Section 4) 
have, effectively, shown that these results are also by-products of the Wiener-Hopf factorization 
met in the theory of random walk, thanks to Kaluza (1928, Proposition 1). 

Obviously, GSR asserts that if W and X are independent random variables (r.v.'s) with W 
nonnegative and X exponential, then WX is i.d., while each of Theorem 5 of Shanbhag & Sreehari 
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(1977) and Theorem VI. 3. 13 of Steutel & van Harn (2004) asserts that the result in question remains 
valid irrespectively of whether or not the r.v. W is nonnegative. Moreover, the discrete version of 
GSR is extended by Theorem VI. 7. 10 of Steutel & van Harn (2004), which, in turn, states, implicitly, 
that if X is an integer-valued r.v. with ¥{X = 0} > and the conditional distributions of X given, 
that X is nonnegative and of —X given that —X is nonnegative, are both completely monotone 
(on {0, 1, ...}), then X is i.d. 

In this article, in Section [21 we present two auxiliary lemmas and, in Section [31 we extend the 
argument based on the Wiener-Hopf factorization, met in C.R. Rao et al. (2009), to establish certain 
results based on log-convex sequences and log-convex functions, respectively, and demonstrate, also, 
that Theorems VI. 7. 10 and VI. 3. 10 of Steutel & van Harn (2004) are their corollaries. Besides, we 
include, in both of these sections, various remarks with comments, comprising examples in some 
cases, on our main findings as well as on results in the existing literature, that are of relevance to 
the present investigation. 



2 Auxiliary lemmas and related observations 

We may start this section by stating first the following crucial definition: 

Definition 2.1 Let B C R such that B = {a, a + 6, a + 26, . . .} or B = (a, oo) with a € K and 
6 £ (0, oo). Then, any function g : B i— > [0, oo), such that 

g 2 (x + y) < g(x) g(x + 2y), for each x£B and y&C, (1) 

where 

c _{{b} if B = {a, a + 6, a + 26,...}, 
\ (0,1) if B = (a,oo), 
will be referred to as a log-convex function (on B). 

(Note that, for notational convenience, in Definition 12.11 we allow g, that is identically equal to 0, 
to be also called log-convex, and that Zygmund (2002, p. 25) has used the same convention in his 
definition of a log-convex function.) 

Obviously, by Definition 12.11 a log-convex function g that is not identically equal to on 
B \ {a} is non- vanishing, and, essentially, by (an observation in) Loeve (1963, p. 159), in the case 
of B = (a,oo), a log-convex function (on B), as defined in Definition 12. 1} is either continuous or 
not a Borel function; for the literature supporting the claim of Loeve (1963, p. 159), one may 
consult, e.g., Hardy et al. (1952, p. 96) or Donoghue (1969, Section 1.3). It may also be pointed 
out that, if g is a log-convex function, then given any x and y as in ([Tjh we get the restriction to 
{x, x+y, x+2y, . . .} as log-convex on that set, and we refer to the sequence {g(x+ny) : n = 0, 1, . . .} 
as a log-convex sequence. 

In what follows, in this section, we give the auxiliary lemmas and remarks of relevance to these. 



3 



Lemma 2.1 Let B be as in Definition \2.1\ let (Q,£,fi) be a measure space, and let {g x : x G B} 
be a family of integrable functions on the measure space so that, for almost all (for short, a. a.) 
uj 6 O, g x (u>), x G B, are log-convex. Define (in the notation of Loeve (1963, p. 119)), 

h(x) = j g x , x G B. 

Then, h is a log- convex function (on B). (We do not impose here the restriction that h be continuous 
if B = (a,oo).) 

Proof. If a, /3, 7 G (0, 00), then, for instance, by the first statement in the last paragraph of 
Kingman (1972, p. 18), we have 

(3 2 < a 7 (2) 

to be equivalent to that 

aX 2 - 2/3A + 7 > 0, for all A > 0. (3) 

In view of the assumptions of the lemma and the criterion of ([!]) for a function to be log-convex, 
it follows that if x G B and y G C (as in dU), we have ([3]) to be valid with a, (3 and 7 replaced, 
respectively, by g x (w), 9x+y{u) and g x +2y(^>) for a.a. to G fl, and, hence, to be valid with a, f3 and 
7 replaced, respectively, by h(x), h{x + y) and h(x + 2y). This, in turn, shows, since ([3|) implies 
([2]), that (H]) holds with h in place of g, and, consequently, that the lemma holds. □ 

Lemma 2.2 Let G be a Lebesgue-Stieltjes measure function (i.e., a nondecreasing, right- 
continuous, real-valued function) on R, so that, for some a G R ; G{x) = if x < a and G(x) 
is differentiable with log-convex derivative if x > a. Define the sequence {G n : n = 1,2,...} of 
Lebesgue-Stieltjes measure functions on R, so that, for each m,n G {1, 2, . . 

( G(a + m/n) if x G [a + (m — l)/n, a + m/n), 
G n (x) = < 

I if x G ( — 00 , a). 

Then, {G n : n = 1, 2, . . .} converges weakly to G and, for each n G {1,2,...}, G n is concentrated 
on B n with G n (x) — G n {x—), x G B n , log-convex, where B n = {a, a + 1/n, a + 2/n, . . .}. 

Proof. That {G n : n = 1, 2, . . .} converges weakly to G is obvious, since, by assumptions, for each 
nG{l,2,...}, 



\G n {x)-G{x)\=G n {x)-G{x) 



< G(x + 1/n) - G(x) if xe[a,oo), 
= if x G (—00, a). 



It is also obvious, by assumptions, that, for each n G {1, 2, . . .}, G n is concentrated on B n and, in 
view of the observation below 11.82 in Titchmarsh (1978, p. 368), that 

/•1/n 

G n (x) - G n (x-) = G(a)S(x) + / G'{x + y)dy if x G B n , (4) 

Jo 
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where S(a) = 1 and S(x) = if x > a. Clearly, we have in flU), the function G(a)S(x), x € B n , to 
be log-convex, and, in view of the log-convexity of G' on (a, oo), for each y € (0, oo), the function 
G'(x + y), x € B n , to be log-convex. Consequently, by Lemma [2TTI it follows that G n {x) — G n (x—), 
x € B n , is log-convex, and, thus, the lemma holds. □ 

Remark 2.1 In the literature, usually, a log-convex sequence {u n }, with uq = 1 and < u n < 1 
if n > 0, is called a Kaluza sequence, especially, in recognition of the findings on such sequences 
given in Kaluza (1928); to be short, we refer to this as a KS. Indeed, Proposition 1 of Kaluza (1928) 
establishes that each KS is renewal. Kingman (1972, Section 1.5), in particular, and Shanbhag 
(1977) amongst others, have made some further observations on these sequences. [Incidentally, the 
statement in Kingman (1972, p. 18) that the class of KS's is a closed convex subset of the class of 
renewal sequences requires the convention of the sequence {u n }, with no = 1 and u n = if n > 0, 
being Kaluza, adopted. A similar blemish, also, exists in the KS-related proof, given by Goldie 
(1967) for his Theorem 2, since it assumes, implicitly (in the notation used in it) that ¥{Z = 0} < 1, 
without stating that this is so.] If g is a non- vanishing log-convex function on B, then ([TJ implies 
that, for each pair x and y as in ([1]), the sequence {g{x + ny)jg(x + (n — l)y) : n = 1, 2, . . .}, is 
nondecreasing, and, hence, for n € {1, 2, . . .}, 

gjx + ny) = -A- g(x + my) < yr g(x + my) ^ g(x + 2ny) 
9(x) JL1 g(x + (m- l)y) ~ m ±^ +l g{x + (m - l)y) g(x + ny) ' 

implying that 

g 2 (x + ny) < g(x) g(x + 2ny). (5) 

In view of ([5]), it follows that ([TJ is equivalent to its version with C replaced by C* , where 
C* = {6,26,...} in the discrete case, and C* = (0,oo) otherwise. Obviously, there are several 
other equivalent formulations for ([TJ) . 

Remark 2.2 In view of the result of Titchmarsh (1978) met in the proof of Lemma [2.21 it follows, 
by Lemma 12. 1\ that {n(G(x + {1/n)) — G(x)),x > a : n = 1,2,...}, with G as in Lemma 12.21 
is a sequence of continuous log-convex functions; since the sequence converges (pointwise) to the 
function G'{x) , x > a, it is seen, by (10.4) and (10.7) in Volume I, Chapter 1, of Zygmund (2002), 
that this latter function is indeed a continuous log-convex function. That G' is continuous, follows, 
also, from Loeve (1963, p. 159), since it is a log-convex Borel function; we could have, obviously, 
used this information, in place of that in Titchmarsh (1978), in the proof of Lemma [2. 2 i It may be 
noted, in this connection, that Loeve (1963) and Zygmund (2002) use different approaches to define 
convex functions, though, these turn out, in view of (10.7) of Zygmund (2002), to be equivalent 
in the case of continuous functions defined on open intervals, and we follow (in Definition 12. 1|) 
implicitly Loeve (1963) to define a log-convex function g in the case of B = (a, oo). 
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Remark 2.3 (i) If B = {a, a + b, a + 2b, . . .} with a£l and b G (0, oo), then, by Lemma [2TTT any 
nonnegative real function h on S, of the form 

oo 

/i(x) = <?( x + n b)v n , x&B, (6) 

n=0 

where g is a log-convex function on Z? and {v n : n = 0, 1, 2, . . .} is a sequence of nonnegative reals, 
is log-convex. 

(ii) Suppose B = (a, oo) with a G R, and 5 is a continuous log-convex function on B. Then, 
in view of Lemma |2.1( by Loeve (1963, p. 159), it follows that (since it is a Borel function) any 
nonnegative real function on B, satisfying 

h(x) = g(x + y)du(y), x G B, (7) 

J[0,oo) 

with v as a Lebesgue-Stieltjes measure, is a continuous log-convex function. Similarly, if a 
nonnegative, using the relevant information (i.e., relative to C*) appearing in Remark 12. 1\ it can be 
seen that any nonnegative real function h on B, satisfying (|7|), with "x + y" replaced by "xy" and 
"[0,oo)" replaced, respectively, by "(0, 00)" if a = and "[1, 00)" otherwise, is also a continuous 
log-convex function. (Since, in each case, {h n : n = 1,2,...}, where, for each n, h n satisfies 
the relevant version of ([7]), with v(- R [0, n]) in place of v, is a sequence of continuous functions 
converging, by the monotone convergence theorem, to h, that h has the stated property is implied, 
in view of Lemma 12.11 and the last observation in Remark 12.21 also, by Zygmund (2002; Theorem 
10.4).) 

(Hi) If a is nonnegative and H is a Lebesgue- Stieltjes measure function on R, such that 

f E(G(xV)) if x > a, 
H(x) = { 1 1 » - (8) 

I otherwise, 

where G is as in Lemma 12.21 and V is a nonnegative real r.v. meeting, if a > 0, a further condition 
of ¥{V > 1} = 1, it follows (on reading, for convenience, 0G'(0) = 0), by Fubini's theorem, that 

H{x)-H(a)= I E(VG'(yV))dy, are (a, 00). (9) 

J (a,x) 

Clearly, then, in view of the properties of G' referred to in Remark 12.21 and our observation on 
the latter version of h in (ii) above, we have that H is differentiable on (a, 00) with continuous 
log-convex derivative (E(VG'(xV)) for x > a). 

Remark 2.4 Proposition 1 of Kaluza (1928) obviously implies that any nondegenerate log-convex 
probability function on {a, a + 6, a + 26, . . .}, with a G R and b G (0, 00), corresponds to an r.v. 
a + bX, where X is an r.v. with compound geometric distribution, and, hence, is i.d. In view of 
this, Lemma [2 . 2 1 implies that if G in it is a d.f. or H of Remark 12.31 (Hi) is a d.f., then it is i.d.; this 
latter result is an extension of Theorem 4.2.6 of Steutel (1970) (or of Theorem III. 8.4 in Steutel 
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and van Harn (2004)). The two results based on Kaluza's proposition, met here, clearly extend 
(respectively, in the reverse order) GSR and its discrete analogue. We have already given some 
information to this effect in the introduction. In view of this, C.R. Rao et al. (2009, Section 4), 
tells us that Theorem VI. 3. 13 and Theorem VI. 7. 10 of Steutel & van Harn (2004) have alternative 
proofs based, at least partly, on KS's or, in particular, completely monotone sequences. 



3 Results based on log-convex sequences and log-convex functions 

To recall partially the information already implied, in our discussion so far, especially in Remark 
12.41 applying the Wiener-Hopf factorization relative to a random walk, C.R. Rao et al. (2009, 
Section 4) effectively proved Theorems VI. 3. 13 and VI. 7. 10 of Steutel & van Harn (2004). As 
hinted in the remark referred to, this approach relies also upon certain properties of log-convex 
sequences, or, in particular, of KS's. Included in the following two subsections, viz., Subsections 
13.11 and 13. 2\ of the present section, are some further KS-related results and remarks. In Subsection 
13.11 modifying the relevant arguments in C.R. Rao et al. (2009) appropriately, we establish a 
key theorem, viz., Theorem 13. 1\ with assertion in it based on log-convex sequences, and, then, 
give four of its important corollaries and, in Subsection 13. 2\ make some pertinent observations 
through remarks on the results obtained. Prom the presentation of Subsection 13.11 it is clear that 
Corollaries 13.11 and 13.31 with assertions based on log-convex sequences and log-convex functions, 
respectively, imply Corollaries 13.21 and 13.41 respectively, and, also, that the latter corollaries are 
respective rephrased versions of the aforementioned Steutel & van Harn (2004) theorems. 

3.1 The key theorem and its corollaries 

We now present in this subsection our main result, viz., Theorem 13.11 and four of its corollaries, 
referred to above. 

Theorem 3.1 Let F be a d.f. relative to a probability distribution {p x : x = 0, ±1, ±2, ...} such 
that, for some constant K > 0, 



Y^jL\ x \+i v -j if x — -1, -2, . . . , 

YSU+iVj if a: = 1,2,..., (10) 

max {E^i^-i>E^i^} if x = °i 



with {vj : j = ±1, ±2, . . .} as a sequence of nonnegative reals for which {«-•/ : j = 1,2, . . .} and 
{vj : j = 1,2,...} are log-convex. Then, F is i.d. 

Proof. We may assume, without loss of generality, that {vj : j = 0, ±1, ±2, . . .}, with vq denoting 
some nonnegative real number, is so that Yl'jLi v -j = Y^ < jL\ v j an£ f Ylj=-i v j = ^ this follows 
since if YljLi v -j YlfLi v j-> then, taking, without loss of generality, Y^jLx v -j > Y^=i v j-> we 
can replace v\ by v% + Y^T=i v -j ~ YlTLi v j an< ^ verify that the normalized version of the resulting 



7 



■u-sequence is as required and (jlOp (possibly, with different K) holds with the original v-sequence 
replaced by this. Now, since we need a proof for the theorem only when the distribution (satisfying 
(jlOp ) is nondegenerate, we can assume, again, without loss of generality, that Vj > for all j > 0, 
and define, for each positive integer k, a sequence {v^ : j = 0, ±1, ±2, . . .} satisfying 



(fe) 



if j = 0, -1, ±2, . . . , ±k, —(k + 1), —(k + 2) 



^(^)' " if J = fc + l,A; + 2,..., 



Kingman (1972, p. 18) involves the idea that we have used in the construction of {uj } from 
{vj} and implies, in view of (fTOj) . that < v^+i/vk < 1; note, in particular, that, in the present 
case, {ui+j/wi : j = 0,1,2,...} and {v-i-j/v-i : j = 0,1,2,...} are such that the first one is 
a decreasing KS and, unless it is the sequence {1,0,0, . . .}, so also is the second one. Obviously, 
for each positive integer k, {{v_l + Vq + u£ ) -1, yj : j = 0, ±1, ±2, . . .} is a sequence possessing 
the properties sought of {vj : j = 0, ±1, ±2, . . .} above and provides us with a distribution {pi^} 
satisfying the version of (jlOp (possibly, with different i^T and) with Vj's replaced by Vj^'s- Denoting 
the d.f. relative to this latter distribution by F^, it is easily seen that {F^ : k = 1, 2, . . .} converges 
weakly to F, and hence, in view of the closure property of the class of i.d. distributions under weak 
convergence, it is sufficient if we prove that, for each positive integer k, F^ is i.d. 

Consequently, it follows that to prove that F is i.d., there is no loss of generality if we 
assume that for some jo € {1,2,...} and b S (0,1), Vj oc b 3 , j > jo. Assume then that 
{vj : j = 0, ±1, ±2, . . .} meets this additional condition and define c = — In b. Also, define now the 
sequence {wj : j = 0, ±1, ±2, . . .} so that 

( v i ~ v j-i if 3 = - 1 ) -2,... , 

Vj-Vj+i if j = l,2, (11) 
vo if j = 0, 

and observe that this is a distribution with mean 0, having an m.g.f. with its domain of definition 
as a superset of [0, c). Obviously, in the present case, {p x } has the same property as {vjj} that it 
has an m.g.f. with [0, c) as a subset of its domain of definition; denote then the m.g.f. 's of {p x } and 
{wj} by M and M*, respectively. Following the relevant approach of C.R. Rao et al. (2009, Section 
4.2) (involving Theorem XII. 2. 2 of Feller (1971) or otherwise), we can now conclude that the weak 
descending and the ascending ladder height measures associated with the random walk relative to 
{vjj} are indeed probability measures with m.g.f. 's having domains of definitions as supersets of 
[0, c), and that 

(1-M*(t)) fl-MUt)\fl-MX(t)\ , . 



;i-e-*)(l-e*) 

where and M| are, respectively, the m.g.f. 's of the weak descending and the ascending ladder 
height measures referred to. (For a simple argument to see that the first equation in (|12p holds, 
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refer to Remark 13.21 (i), given below.) In view of (XII. 3. 7a) in Feller (1971) corresponding 
to the distribution of a weak descending ladder height and its analogue corresponding to an 
ascending ladder height, it follows that, in our case, the distributions referred to are of the 
forms {Y^=o w i-n v n^ ■ i = 0,-1,-2,...} and {£)£L w i+n v n '■ i = 1,2,...}, with z4 and 

(2) 

Vn nonnegative for all n, respectively. Since for any nonnegative integer-valued random variable 
Z, we have a standard result that, for each t ^ for which E(e ) < oo, 

oo 

(l-E(e tz ))/(l-e t ) = Y,e tj nZ>j}, 

on appealing to (|12p . especially, in view of Fubini's theorem, Remark 12.31 (i) and the log-convexity 
properties of {v-j : j = 1,2,...} and {vj : j = 1,2,...}, it is clear that there exist sequences 
{v*j '■ J ' = 0, 1,2, ...} and {v 2 j ■ j = 0, 1, 2, . . .} of nonnegative reals and of positive reals, 
respectively, that are log-convex such that 

KM(t) = (^£ y fj v* 2j e^ , t E [0, c). (13) 

(Note that we have assumed in this proof, without loss of generality, that Vj > for all j > 0.) We 
have, obviously, {u^ - : j = 0, 1, 2, . . .} and {v^j : j = 0, 1, 2, . . .} in (fTHj) to be, in view of the related 
information in Remark 12.11 proportional to renewal sequences and, hence, to appropriate discrete 
i.d. distributions (with the first one as degenerate at the origin or compound geometric, and the 
second one as compound geometric). We have then that F is i.d., implying that the theorem holds. 
□ 

Corollary 3.1 For each a G [0, 1], let F^ denote a d.f. so that 

F (q) (x) = aF 1 (x) + {l-a)F 2 (x), iGl, (14) 

where F± is the d.f. relative to the degenerate distribution at the origin and F 2 is the d.f. of a 
nondegenerate distribution satisfying i!0\). Then, each F^ is i.d. 

Proof. We need a proof only when a E (0, 1). Note now that, in this case, for each a, the 
probability function (on {0, ±1, ±2, . . .}) relative to iA°0 satisfies ()10p . with K unaltered and {vj} 
replaced by {v^}, where 



( a ) _ \ Ka+ {I- a)Vj if j = -1,1, 
[ (1 — a)vj otherwise. 

Consequently, we have the corollary. □ 

Corollary 3.2 Each probability distribution {p x : x = 0, ±1, ±2, . . .}, with pq > and {p~ x '■ x = 
0, 1, 2, . . .} and {p x : x = 0,1,2, .. .} completely monotone, is i.d. 
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Proof. In view of Hausdorff's theorem, referred to, e.g., in Theorem VII. 3.1 of Feller (1971), it 
follows that 



Pi 



PomPj if a? = -1,-2, 

Pomi if x = 1,2,... 

Ai) f _ „ (2h 



{ Po< J (=Po<0 if x = 0, 

where {rn£^ : x = 0,1,2,...} and {m^ : x = 0,1,2,...} (with, obviously, JTtg = = 1) 

are moment sequences relative to probability distributions concentrated on [0, 1). Defining 
{vj : j = ±1, ±2, . . .} for which 



Pj+2 ~ Pj+i if 3 = ~ 2 ' ~ 3 ' ~ 4 > ■ • • j 
Pj-Pi+i if i = 0,1,2,..., 



it is then seen that {p x } satisfies (jlOl) with vj's meeting the required condition and K = 1. Hence, 
we have the corollary. □ 

Corollary 3.3 For each a G [0, 1], Zei 6e so i/tai is met but for that F2, now, instead of 
that in the statement of Corollary \3.1\ is the d.f. relative to an absolutely continuous distribution 
with density f 2 satisfying 

, . I [1°°! v\{y)dy if x < 0, , , 

f 2 (x) = ) J N uy; y ' (15) 

[ f™v 2 (y)dy if x > 0, 
with v\ and V2 as (nonnegative) log- convex functions on (0, 00). Then, each is i.d. 

Proof. For each n E {1, 2, . . .}, define the d.f. F^ such that it is concentrated on {0, ±— , . . .}, 
having {F 2 ((x + l)/n) — (x/n) : x = 0, ±1, ±2, . . .} to be of the form of {p x } of Theorem 13. 11 
satisfying ([10]) with iT n > (where K n — > 1 as n — > 00) in place of K > and {u,} replaced by 
{v^}, where 

v (n) = f Jo /n Io /n M(\j + M/n) + y + z)dydz if j = -1, -2 

'' I Jo /n fo /n V2(((j-l)/n)+y + z)dydz if j = 1, 2, . . . . 



Note that, for each n, given y, 2, G 0, 1/n), |u r ((m/n) + y + z) : m = 0, 1, . . .}, r = 1, 2, are log- 
convex sequences. Hence, by Lemma 12.11 it follows that, for each n, {v^} meets the requirements 

(n) 

of the theorem. Since F 2 can easily be seen to be so that it converges to F% weakly, essentially, 
in view of Corollary 13.11 and the closure property of the class of i.d. distributions, we can conclude 
that, for each a G [0, 1], F^ a ' is i.d., and, thus, we have the corollary. □ 

Corollary 3.4 For each a G [0, 1], let F^ a ' be so that |i^[ ) is met but for that F2, now, instead of 
that in the statement of Corollaru \3.1i is the d.f relative to an absolutely continuous distribution 
with density f% for which, for each r G {1,2}, the function /2((— l) r t), t G (0, 00), is completely 
monotone. Then, each F^ is i.d. 
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Proof. In view of the conditions to be met by /2, essentially, appealing to a version of Bernstein's 
theorem appearing as Theorem XIII. 4. la of Feller (1971), it follows that, in this case, (|14h holds 
for each a G [0, 1], with u r 's in (fT5j) so that 

v r (y) = [ \ 2 e- Xy dn r (\), y G (0,00), r = 1,2, 

■/(O.oo) 

where /i r , r = 1,2, are measures so that /ii + /X2 is a probability measure concentrated on 
(0, oo). Clearly, and V2 in this case are log-convex on (0, oo), satisfying (in obvious notation) 
v r (.)v r (.) — (v r ) 2 (.) > 0, r = 1, 2. Hence, we have the corollary from Corollary 13.31 □ 

3.2 Some relevant remarks 

We devote the present subsection, as implied before, to making specific observations on our findings 
through remarks: 

Remark 3.1 The equation (|10p in the statement of Theorem 13.11 can be also expressed, with a 
minor notational adjustment, so as to have K = 1. However, to have the options such as that 
in which {v-i-j : j = 0, 1, 2, . . .} and {v\+j : j = 0, 1, 2, . . .} are moment sequences open to us, 
and to avoid unnecessary notational complications in the arguments used to prove Theorem 13.11 
we have decided to retain K in the equation referred to. Also, in view, especially, of the relevant 
information in Remarks l2.3l (i). 12.31 (ii) and !2.4l it is obvious that if Vj = for some (and hence all) 
j > 1 or j < —1, Theorem 13.11 and Corollary 13. 11 and, if g\{y) or 52(2/) = for some (and hence all) 
y > 0, Corollary 13.31 follow as simple consequences of Proposition 1 of Kaluza (1928). To illustrate 
that Corollaries 13.21 and 13.41 are more restrictive versions of Corollaries 13.11 and 13.31 respectively, we 
now give the following simple example: 

Example 3.1 Let g be a function defined on (0, 00), such that 

g{x) = e -*+ft(*) j x > 0, (17) 

where h(x) = (1 — x) 2 if x G (0, 1) and h(x) = if x > 1. Then, g and g(y)dy, x > 0, are log- 
convex but not completely monotone; the log-convexity of the first function is obvious and of the 
second function follows by the relevant information in Remark 12.31 (ii). and that the functions are 
not completely monotone follows, since these are not differentiable twice and thrice, respectively, at 
the point x = 1. Also, if we now define {v n : n = 0, 1, 2, . . .} such that vq = e and, for each n > 0, 
v n = g(n/2), where g is as in (|17p . then, by the Hausdorff theorem, neither of the sequences {v n } 
and {Ylm=n+i v m-i '■ n = 0, 1, 2 . . .} is completely monotone; note that the cited theorem implies 
that a real sequence {u* : n = 0, 1,2, . . .}, for which eu* +2 — 2e 1//2 t>* +1 + = for some n > 1, 
turns out to be completely monotone, only if e v% — 2e 1 / 2 ^2 + v* = (and this criterion is not met 
for the two sequences). However, it is easily seen now that these sequences are indeed log-convex. 
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Remark 3.2 (i) A simple argument to show that the first equation in (|12p is valid, is as follows: 
Clearly, in view of the assumptions in the statement of Theorem 13. II in conjunction with those that 
we added in the proof of the theorem claiming that there was no loss of generality in doing so we 
have from (jlOp (in obvious notation) 

#"(l-e-*)(l-e*)M(t) = K(2-e t -e~ t )M(t) 

oo 

= e tx [K(2p x -p x -i -px+i)] 



-oo 

oo 



1- e tx w x = (l-M*(t)), te[0,c). 



-oo 



K(2p x -p x -i -Px+i) 



This is obvious on noting, in particular, that the last but one equation holds since 

v x -i-v x = -w x if x = -1,-2,..., 

v x+ i-v x = -w x if x = l,2,..., 

V-\ + v\ = 1 — vo = 1 — wo if x = 0. 

(ii) The argument that we have used in the proof of Theorem 13.11 remains valid with vq = 0. 
However, to make the link between our approach in this case and that appearing in Section 4.2 of 
C.R. Rao et al. (2009) more transparent, we have allowed here also the case vq ^ 0. To illustrate 
this, we now consider the following example. 

Example 3.2 Let n € {1,2,...,} and {p x : x = 0,±1,±2, ...} be a probability distribution, 
discussed implicitly in Section 4.2 of C.R. Rao et al. (2009), so that 

Es=i c isPiT if x = -1,-2,..., 
1>< - S E n s =iC2sP x 2s if x = l,2,..., 
max{ c\,c^} if x = 0, 

where c* = J2™=i c rs, r = 1,2, c rs > and p rs G (0, 1) for each r G {1, 2} and s 6 {1,2,..., n}. 
Then, by Corollary 13.21 it is immediate that {p x } is i.d. (Note that, in this case, the m.g.f. for 
the distribution exists, with domain of definition having as an interior point, and hence there are 
obvious advantages.) However, we may stress here that the specific construction that C.R. Rao et 
al. (2009) have given to get M* from M assumes (in our notation) vo = Wq > (as implied by 
(4.4) of the cited paper). 

Remark 3.3 From the proofs that we have given above, especially for Theorem 13. II and Corollary 
13.31 it is obvious that, for each of Theorem 13.11 and Corollaries I3.HI3.41 there exists a sequence 
{Xi^ n — X2 t n ■ n = 1,2, . . .}, with, for each n, Xi >n and X2^ n as independent discrete nonnegative 
r.v.'s having log-convex distributions (on {0, b n , 2b n , . . .}, for some b n > 0), converging in 
distribution to an r.v., with distribution, that is claimed to be i.d., in the respective assertion; 
from C.R. Rao et al. (2009, Section 4), it follows further that, in the case of Corollaries 13.21 and 13.41 
the observation remains valid with "log-convex distributions" replaced by "completely monotone 
distributions or their scale variations". 
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Remark 3.4 It follows easily that the classes of the mixtures met in the statements of Corollaries 
13-H 13.31 and 13.41 and the class of the distributions in the statement of Corollary 13.21 are all convex 
(with members that are i.d). In view of this, it is obvious, amongst other things, that if X is an 
r.v. with distribution as a member of the class of mixtures in the statement of Corollary 13.31 (or, 
in particular, in the statement of Corollary 13. 4p and W is an r.v. independent of X, then, by the 
closure property of the class of i.d. distributions, WX and W|X| are i.d. (Moreover, in view of 
the relevant information in Remark 12.31 (ii) or in the proof of Corollary 13.41 respectively, Fubini's 
theorem implies that the distributions of WX and W|X| lie in the class relative to the distribution 
of X, referred to.) It may be worth pointing out in this place that if X* is an r.v. with its d.f. as 
G of Remark 12.41 but with o = 0, and W is a nonnegative r.v. independent of X*, then, by (the 
related information in) Remarks 12.31 (in) and 12.41 WX* is i.d, since, there is no loss of generality 
in assuming W > and we have the d.f. of WX*, in this case, to be a specialized version of H of 
Remarks (in), with V = 1/W. 

Remark 3.5 In the previous remark, we came across some cases of independent r.v.'s W and X 
with WX i.d. One may now raise a question as to whether discrete versions of these relative to 
integer- valued r.v.'s are valid. That the answer to this question is in the negative is shown by the 
following example. 

Example 3.3 Let W and X be independent r.v.'s with X geometric and W a nondegenerate r.v. 
whose distribution is concentrated on {2,3}. Then, since WX is a nonnegative integer-valued 
r.v. with support of its distribution so that it includes the points 0, 2 and 3, but not the point 
5, Theorem 4.2.3 of Steutel (1970) or Theorem II. 8. 2 of Steutel & van Harn (2004) implies that 
this latter r.v. is non-i.d. (In view of the closure property of the class of i.d. distributions, one 
can, obviously, obtain more general examples to illustrate that this is so.) Incidentally, the present 
example illustrates also that there exist non-i.d. mixtures of compound geometric distributions. 

Remark 3.6 Let X and X* be r.v.'s as in Remark I3T41 Also, let (f> : K >->• M and (f>* : [0, oo) h-> [0, oo) 
be one-to-one and onto (i.e., bijective) functions with 0(0) = <^*(0) = 0, for which their inverses 
and (<fi*)~ l are such that 

= | -Jt l Mv)dy if x < o, 

1 Joip2(y)dy if x > 0, 

where ip r , r = 1,2, are decreasing log-convex functions on (0, oo) with (in standard notation for 
derivatives) —ip r (y), y € (0, oo), also as log-convex for r = 1,2, and 

(cf ) *)- 1 (x)= f r(y)dy if ^>0, 
J o 

where ifi* is a decreasing log-convex function on (0, oo). Then, denoting respectively by H and H* 
the d.f.s of 4>(X) and <f)*(X*), it is seen, for example, that, for x > 0, (in standard notation for 
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derivatives) —H"(x) and {H*)'(x) are both (nonincreasing) continuous log-convex functions. (To 
understand this last claim properly, note that if g is proportional to a log-convex density on (0, oo), 
then, by Remark 12.31 (n), g{{<j))~ l (x)) , x > 0, g(((p*)~ 1 (x)), x > 0, and g(y + 4>~ 1 (x)) dy, x > 0, 
are (nonincreasing) continuous log-convex functions, and, also, that the products (by Definition 
12. 1[) and the sums (by Lemma 12. 1 j) of finitely many log-convex functions on B are log-convex on 
B.) Hence, in view of what we have already pointed out in Remark 13.44 it follows in particular 
(appealing to the underlying symmetry in the case of H) that if W\ and W2, with W2 nonnegative, 
are r.v.'s, independent of X and X* , respectively, then, for all n € {1,2,...}, W\X n and, for all 
f3 € [l,oo), and W%{X*)^ are i.d. (Clearly, there exist several other possibilities here.) 

Remark 3.7 If {p x : x = 0, 1, . . .} is a probability distribution so that p x = u x — u x+ i, x = 0, 1, . . ., 
with {u x : x = 0, 1, . . .} completely monotone, then, by Corollary 13.21 it is i.d. However, that this 
assertion does not remain valid, if we replace the condition that {u x } is completely monotone by 
that it is a KS, is shown by the following example. 

Example 3.4 Let 6, c £ (0, 1) with b < c, and {u x : x = 0, 1, . . .} be so that 

f b x if a; = 0,1, 

{ b 2 c"- 2 if £ = 2,3,... . 

Also, let p x = u x — u x+ i, x = 0, 1, . . ., and / be the corresponding probability generating function 
(p.g.f.) . Observe now that {u x } considered is a Kaluza sequence and if we choose c sufficiently 
close to 1 (with b fixed), then, {p x }, in this case, turns out to be non-i.d., since, in spite of that 
< f(0) < 1, we cannot have here /(•) = e~ x+Xg ^ with A > and g as a p.g.f. satisfying g(0) = 0. 
(To see the validity of the claim on {p x }, it is sufficient if we verify that ln/(s) has its second 
derivative at s = to be negative.) 

Remark 3.8 One may now raise a question as to whether Theorem 13.11 remains valid if we take 
in place of ()10p its version, obtained from it, replacing the four summations, on its right hand side, 
respectively, by v x +i, ^-i an d v%. Taking a hint from Example 13.41 we can now construct 

the following example to show that the answer to the question is in the negative. 

Example 3.5 Let {p x : x = 0, ±1, ±2, . . .} be a probability distribution satisfying 

( c l*l if x = 0,-1,-2,..., 

b 2 c x ~ 2 if x = 2,3,..., (18) 
b if x = 1, 



K Px 



with K > and b, c £ (0, 1) so that b < c. Clearly, (|18p may, now, be viewed as the version of (|10p . 
that is sought, with, e.g., v\ = v-\ = 1 and Vj's for j ^ —1, 1 defined in obvious way. Consequently, 
if X is an r.v. whose distribution is {p x } then (with t complex) 

KE(e tx ) = (l-ce-y' + be' + ^e^il-cey 1 , 

= (1 - 6c) (1 - ce-*)"Ve*), Re(t) <E (lnc, — lnc), (19) 



14 



with g(s) = 1 + bs + b 2 s 2 (l — c 2 )(l — 6c) -1 (1 — cs) _1 , for each s = e*, Re(i) G (lnc, — lnc). Suppose 
now that we specialize to the case with b + b 2 (l + c)(l — fcc) -1 < 1 and (1 — c 2 )(l — 6c) _1 < 1/2 
(which subsumes, e.g., the case with b = 1/3 and c = 8/9). On applying the standard power-series 
expansion for ln(l + z), \z\ < 1, by minor manipulation, f|19|) implies then that, in this special case, 
there exists a sequence {q x : x = 0, ±1, ±2, . . .} of reals, with q$ = 0, ^2 < 0, = c'^'/W if ^ < 0, 
and Y^=2 I I < 00 > sucri that 

oo 

ln(E(e tx )) = £ (e te - l)q x , Re(i) G (lnc,0], 

£ = — OO 

asserting, in view of Lukacs (1970, Remark 1, p. 118) or Blum & Rosenblatt (1959, Theorem 1), that 
{p x } is non-i.d. (Incidentally, Blum & Rosenblatt (1959, Theorem 1) tells us that any nondegenerate 
i.d. distribution with at least one discontinuity is so that its characteristic function (ch.f.) has a 
unique representation as the product of a degenerate ch.f. and a ch.f. that is of the form e~ x+Xip ^ 
with A > and ip as the ch.f. of a d.f. G satisfying G(0) — G(0— ) = 0; note that if the i.d. 
distribution corresponds to an integer- valued r.v., then the degenerate ch.f. and G referred to here 
also correspond to such r.v.'s.) 

Remark 3.9 A further example, i.e., Example I3.6| that appears below, tells us that Corollary 
13.31 does not remain valid if "Jj^? v\(y)dy n and U J^° V2{y)dy v in (|15p are replaced by "«i(|x|)" and 
u V2(x)" , respectively. (Incidentally, as a by-product of this, the argument used in the proof of 
Corollary I3.3[ with obvious alterations, then implies that the answer to the question in Remark l3.8l 
is in the negative, thus, supporting the conclusion of Example 13.51 ) 



Example 3.6 Let F be an absolutely continuous d.f. with p. d.f. / such that, for an appropriate 
constant K > 0, 

, f e- 5 ^ if x < 0, 

K f( x ) = \ - Sx+h(x) . f ~ ( 20 ) 
[ ae dx + n ( x > if x > 0, 

with < a < e^ 1 , S > and h as in Example 13.11 If X is an r.v. with d.f. F, we see then, using, 
in particular, Fubini's theorem or the method of integration by parts, that (with t complex) 

KE(e tx ) = (d + ty^l + aid + t) J™ e tx - 5x+h ( x Ux} 

= (5 + i) _1 {l -ae + a J°° e (*+*)*-a&+fc(*) (25 - h'(x))dx} 

= {l-ae)(5 + t)- 1 g*(t), Re(i) G (—5, 5), (21) 



where 



/•oo 

g *(t) = 1 + a(l - ae)' 1 / e **-&+fc(*) (2<f _ h'(x))dx, Re(t) G (-5,5). 

Jo 



Suppose we now take, for convenience, a(l — ae) 1 < 1/4 and 5 = a 2 , and assume that F is i.d. 
with Levy measure v. Clearly, in this case, < g*(0) — 1 < a(l — ae)~ 1 (5(e — 1) + e + 1) < 1, 
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and the function 0, defined by (f>(s) = (g*(is) — l)/(g*(0) — 1), s £ R, is a ch.f.; denote by G the 
d.f. relative to <f>. (The inequalities for g*(0) — 1, an integral over (0, oo), are obvious, in view of 
the assumptions for a and 5, on expressing it as the sum of the appropriate integrals over (0, 1) 
and [l,oo), respectively.) The table on spectral measures for certain i.d. distributions given by 
Lukacs (1970, p. 120), and the power-series expansion for ln(l + z) with z = (g*(0) — l)cp(s), s € R, 
in conjunction with the uniqueness theorem for the Fourier transform of a finite signed measure 
(given, e.g., as Corollary 1.1.2 in Rao & Shanbhag (1994, p. 2)), imply then, by (|2"Tj) . that 

u((-oo,x))= I \y\~ l e Sy dy, if xG(-oo,0), 

J — oo 

and 

oo 

u((0,x}) = J^n-H-l)"- V(0) -l) n G n *(x), if xe(0,oo), (22) 

n=l 

where, for each n € {1, 2, . . .}, G n * is the n-fold convolution of G with itself. From (|22p . we see that 
there exists a function o : (0, oo) i— > R with lim y ^o + (o(y)/y) = 0, so that, for sufficiently small a, 

u((l, 2]) = (1/2) a 2 f (e^M - l)e h ^h'(y)dy + o(a 2 ) < 0, 
Jo 

leading us to a contradiction, and, hence, supporting the claim made in Remark 13.91 

Remark 3.10 To shed further light on the conclusions of Examples 13.51 and 13.61 we may give 
the following relevant information: Extending g and g* appearing in these examples appropriately 
with notation, for convenience, for the extensions respectively as g c and gt, so that their domains of 
definition are the sets of complex numbers, respectively, with moduli lying in (0, 1/c) and with real 
parts lying in (— oo, S), it is seen that, for some (real) t* G (— oo,0), g c (exp{i,s + t*})/g c (exp{t*}), 
sGi, and gt(is + t*)/g : $(t*), s £ R, are ch.f.'s. (The extensions referred to here can be assumed to 
be analytic continuations of their original versions.) If we now allow c and 5 to vary as (distinct) 
members of a c-sequence tending to 1 and a <5-sequence tending to 0, respectively, then the resulting 
sequences of ch.f.'s converge to the ch.f.'s of certain nondegenerate bounded r.v.'s. Clearly, the 
limiting distributions in the two cases referred to are non-i.d., explaining indirectly, as to why we 
have the contradictions in the two examples. 
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